Quasistatically varying log-normal distribution in the 
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Abstract 

Employing data on the assessed value of land in 1974-2007 Japan, we exhibit a qua- 
sistatically varying log-normal distribution in the middle scale region. In the derivation, a 
Non-Gibrat's law under the detailed quasi-balance is adopted together with two approxima- 
tions. The resultant distribution is power-law with the varying exponent in the large scale 
region and the quasistatic log-normal distribution with the varying standard deviation in the 
middle scale region. In the distribution, not only the change of the exponent but also the 
change of the standard deviation depends on the parameter of the detailed quasi-balance. 
These results are consistently confirmed by the empirical data. 
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1 Introduction 

. Log-normal distributions are frequently observed not only in natural phenomena but also 

■ in social ones. For representative example, the probability density function P{x) of personal 



income or firm size x is considered to obey the log-normal distribution [l]-[2] 

In^ (x/x) 



C3 ■ Pln(.x) = — -== exp 



2^2 



(1) 



in the middle scale region.^ Here ic is a mean value and is a variance. A large number 
of persons or firms are included in the middle scale region. The study of the distributions is 
significant. 

The simplest model which describes the log-normal distribution is the pure multiplicative 
stochastic process defined by 

x{t + 1) = R{t) x{t) , (2) 



*e-mail address: ishikawa@kanazawa-gu.ac.jp 

^In the large scale region, the distribution of personal income or firm size follows power-law [3]. We discuss 
this in the next section. In the low scale region, several distributions are proposed (see Refs. [4]-[5] for instance). 
We do not discuss them in this study. 
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where R{t) is a positive random variable. By applying this process iteratively, we obtain 



x{t) = R{t -1) R{t-2) ■■■ R{0) x{0) . (3) 

The logarithm of this equation is 

log x{t) = log R{t - 1) + log - 2) + • • • + log R{0) + log x(0) . (4) 

If log x(0) is negligible compared to log x{t) in the limit t ^ oo and log R{i) {i = 0,1, - ■ ■ ,t — 1) 
are independent probability variables, logx{t) follows the normal distribution in the limit. As a 
result, this model leads the stationary log-normal distribution Pln{x)- 

Eq. (2) means that the distribution of the growth rate R{t) does not depend on x{t). This 
is known as Gibrat's law [1] that the conditional probability density function Q{R\xi) of the 
growth rate is independent of the initial value xi: 

Q{R\xi) = Q{R) . (5) 

Here x\ and X2 are two successive incomes, assets, sales, profits, the number of employees and 
so forth. The growth rate R is defined as i? = X2lx\ and Q{R\xi) as 

Q(RM = ^5^^ (6) 

by using the probability density function P{xi) and the joint probability density function 
Pir{xi,R). 

As far as firm sizes in the middle scale region, however, it is reported that the growth rate 
distributions do not follow the Gibrat's law (5) (see Refs. [6]-[8] for instance). The log-normal 
distribution in the middle scale region cannot be explained by the pure multiplicative stochastic 
process model (2). Instead, we have shown that the log-normal distribution can be derived [9] 
by using no model such as the pure multiplicative stochastic process. In the derivation, two laws 
are employed which are observed in profits data of Japanese firms. One is the law of detailed 
balance which represents symmetry in a stable economy [11]. The other is a Non-Gibrat's law 
which describes a statistical dependence in the growth rate of the past value [9] . 

In Ref. [10], the Non-Gibrat's law and the static log-normal distribution in the middle scale 
region are uniquely derived from the detailed balance. The derivation has been confirmed by the 
empirical data analysis. In this study, we extend the derivation by replacing the detailed balance 
with the detailed quasi-balance proposed in Ref. [12] to derive a log-normal distribution in the 
quasistatic system. By this procedure, the log-normal distribution is described as quasistatic. 
The derivation is consistently confirmed employing data on the assessed value of land in 1974- 
2007 Japan. 
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2 Static log-normal distribution under the detailed balance 

In this section, we briefly review the study in Ref. [10]. As an equihbrium system, we 
investigate profits data of Japanese firms in 2003, 2004 and 2005 which are available on the 
database "CD Eyes 50" published in 2005 and 2006 by TOKYO SHOKO RESEARCH, LTD. [13]. 

Figure 1 shows the joint probability density function Pi2(a;i,X2) of all firms in the database, 
the profits of which in 2004 (xi) and 2005 {X2) exceeded 0, xi > and X2 > 0. The number 
of firms is "232,497". Prom Fig. 1, we approximately confirm the detailed balance which is 
time-reversal symmetry (xi ^ X2) of Pi2{xi,X2) [H]:^ 

Pl2{xi,X2) =Pl2{x2,Xl) . (7) 

Figure 2 shows probability density functions of profits in the database. The distributions 
are almost stable and the following power-law is observed in the large scale region 

P{x) = Cx-"-^ for x>xo , (8) 

where xq is a certain threshold. This power-law is called Pareto's law [3] and the exponent /i is 
named Pareto index. Notice that the Pareto's law does not hold below the threshold xq. The 
purpose of this section is to exhibit the distribution in the middle scale region under the detailed 
balance (7). 

In order to identify a statistical dependence in the growth rate of the past value, wc examine 
the probability density function of the profits growth rate in 2004-2005. In the database, 
we divide the range of xi into logarithmically equal bins as xi G 4 x [iQ^^*^'^^""^^ , lO^"'''''^"'] 
thousand yen with n = 1, 2, • • • ,20. Figure 3 shows the conditional probability density functions 
for r = log]^o ^ the case of n = 1, • • • , 5, n = 6, • • • , 10, n = 11, • • • , 15 and n = 16, • • • , 20. The 
number of firms in Fig. 3 is "20,669", "89,064", "88,651" and "26,887", respectively. 

From Fig. 3, we approximate log^g q{'''\xi) by linear functions of r: 

logio9(^ki) = c{xi) - t+{xi) r for r>0, (9) 
\ogiQq{r\xi) = c{xi) + t-{xi) r for r<0. (10) 

Here q{r\xi) is the conditional probability density function for r, which is related to that for R 
by logio Q{R\xi) = logio qi^lxi) — r — logio(lii 10)- These approximations (9)-(10) are expressed 
as tent-shaped exponential forms as follows: 

Q{R\xi) = d{xi) for R>1 , (11) 

QiR\xi) = d{xi) for R<1 , (12) 

where d{xi) = lO'^'-^^^/ln 10. Furthermore, Fig. 3 shows that the dependence of c(,ti) on xi is 
negligible for n = 9, • • ■ , 20. We assess the validity of these approximations against the results. 

^Similarly, we approximately confirm tlie detailed balance in the joint probability density functions of the 
profits in 2003-2004 and 2003-2005. 
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Figure 4 represents the dependence of t±{xi) on the lower bound of eaeh bin xi = 4 x 
-|^Qi+0.2(7t-i) ^ Pqj, ^ _ 17^ . . . ^ 20, t±{xi) hardly responds to Xi. This means that the Gibrat's 
law (5) holds only in the large seale region of profits [x > xq)."^ In eontrast, linearly 
increases and t^(xi) linearly decreases symmetrically with logj^g^i — 9, • ' ' > 13. From 

Fig. 4, the slops are described as [9] 

X\ 

t±(xi) = t±(xo) :iz a In — . (13) 

Xo 

The parameters are estimated as follows: a ~ for xi > xq, a ~ 0.14 for Xmin < xi < xq,'^ 
Xo = 4 X 10i+°-2(i7-i) ~ 63,000 thousand yen and Xmin = 4 x 10^+°-2(9-i) ~ 1,600 thousand 
yen. Notice that approximations (9)-(10) uniquely fix the expression of t±{xi) under the detailed 
balance [10]. This derivation is included in the proof in the next section. We call Eqs. (11)-(13) 
a Non-Gibrat's law. 

For n = 9, • • • , 20, the dependence of d{xi) on xi is negligible. In this case, the Non-Gibrat's 
law determines the probability density function of profits as follows: 

P{xi) = 0x1-^"+^^ e""^""' ^ for xi > x^in , (14) 

where t+(xo) — t-(xo) ~ /x [14]. This is the power-law in the large scale region (xi > xq) and 
the log-normal distribution in the middle scale one (xmin < a;i < xq). The relations between 
parameters cr^, x in Eq. (1) and a, n, xq are given by a = /Lt = ^ In Figure 5 shows 
that the distribution (14) fits with the empirical data consistently. Notice that the distribution 
cannot fit with the empirical data, if a is different from the value estimated in Fig. 4 (a = 0.10 
or q; = 0.20 for instance).^ 

3 Quasistatic log-normal distribution under the detailed quasi- 
balance 

In the previous section, the log-normal distribution in the middle scale region is exhibited by 
a Non-Gibrat's law under the detailed balance. The resultant profits distribution is empirically 
confirmed in data analyses of Japanese firms in 2003-2004, 2004-2005 and 2003-2005. The 
profits distribution (14) is static, because the derivation is based on the detailed balance (7) 
which is static time-reversal symmetry observed in the system (Fig. 1). 

^Pujiwara et al. [11] prove that the Pareto's law (8) is derived from the Gibrat's law (5) valid only in the large 

scale region and the detailed balance (7). In the derivation, linear approximations (9)-(10) need not be assumed. 
*Here, a constant parameter a takes different values in two regions. This is not an exact procedure. However, 

for firms which are in the large scale region in both years {xi > xo and X2 > xo) or in the middle scale one 
(smin < XI < Xo and Xmin < X2 < xo) , this procedure is exact. In the database, most firms stay in the same 
region. This parameterization is, therefore, approximately valid for describing the probability density function 
P{xi). This is confirmed in Fig. 5. 

^These empirical data analyses are not restricted in the single term 2004-2005. Similar data analyses are 
checked in 2003-2004 or 2003-2005. 
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In the case that an economy is not stable, the detailed balance should be extended to describe 
the state. In Ref. [12], we have derived Pareto's law with annually varying Pareto index under 
the detailed quasi-balance: 

,ax/). (15) 

It is assumed that, in an ideal quasistatic system, the joint probability density function has 
"a xi^ <-> X2" symmetry where ^ is a slope of a regression line: 

logio X2 = logio a^i + logio o ■ (16) 

The detailed balance (7) has the special symmetry 9 = a = 1. Because the detailed quasi-balance 
(15) is imposed on the system, 9 is related to a as follows: 

= l-^logioa. (17) 

Here 10^ is a sufficient large value compared to the upper bound in which 9 and a are estimated. 

In Ref. [12], these results have been empirically confirmed by employing data on the assessed 
value of land [15]-[16] in 1983-2005. In the derivation, under the detailed quasi-balance, we 
have used the Gibrat's law (5) valid only in the large scale region without linear approximations 
(9)~(10). The purpose of this section is to show that the approximations uniquely fix a Non- 
Gibrat's law even under the detailed quasi-balance (15). After that, we identify the quasistatic 
distribution not only in the large scale region but also in the middle scale one. 

By using the relation Pi2{xi,X2)dxidx2 = Pir{xi, R)dxidR under the change of variables 
{xi,X2) ^ {xi,R), these two joint probability density functions are related to each other 

Pir{xi,R) = a Pi2{xi,X2) , (18) 

where we use a modified growth rate R = X2/{a xi^) . Prom this relation, the detailed quasi- 
balance (15) is rewritten as 

P,Rix,,R) = R-'PiRi(^^y\R-') . (19) 

Substituting Pir{xi,R) for Q{R\xi) defined in Eq. (6), the detailed quasi-balance is reduced to 
be 

Pjxi) ^ l Q{R-^\{x2/a)'/') 

P{{x2/af") R Q{R\xi) • ^ ^ 

In the quasistatic system, we also assume that Q{R\xi) follows the tent-shaped exponential 
forms (11)-(12).^ Under the approximations, the detailed quasi-balance is expressed as 



P{{x2/af/') 



^t+{xi)-t-{{x2/aY"')+l ^21) 



^For instance, the tent-shaped exponential forms are observed in data on the assessed value of land in Japan [15] . 
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for R > 1. Here we use the notation P{x) = P{x)d{x). By setting R = 1 after differentiating 
Eq. (21) with respect to R, the following differential equation is obtained: 



e 



t+{x) - t-{x) + ll P{x) + xP '{x) = , 



(22) 



where x denotes xi. The same equation is obtained for R < 1. 

Similarly, from the second and third derivatives of Eq. (21), the following differential equa- 
tions are obtained: 



t+ (x) + 1- (x) = , t+ (x) + x 



. 



(23) 



The solutions are uniquely fixed as 



t±{x) = t±(xo) ± a In ■ 



XQ 



(24) 



This is the same expression under the detailed balance. 

By using the Non-Gibrat's law (11)-(12), (24) and the differential equation (22), probability 
density functions Pi{xi), P2{x2) are also uniquely reduced to 



Pi{xi) 

P2{X2) 



with 



Ci xi-f"^-^ exp 
Ca X2-^^-^ exp 



1 2 

am — 

xo. 



-e aln^ 



(x2/a 



H2 + 1 

in the large or middle scale region where the dependence of d{xi) on xi is negligible. 
Here we consider two log-normal distributions in the middle scale region: 



(25) 
(26) 

(27) 



^tNi(xi) -- 



1 



exp 



xi\/27rai^ 
1 

■ , exp 



In^ (xi/xi) 
2^? 

In^ (X2/X2) 

2(722 



(28) 
(29) 



;/i — ^ In ^ - 

' ^1 ~ ctP Si ' ft 



XI ' 



By comparing Eqs. (25)-(26) to (28)-(29), we identify 0a 
II2 = In • Consequently, the relation between a\ , 02 and Q is expressed as 



and 

2(T2 



(30) 



This is the equation which quasistatic log-normal distributions satisfy. 
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4 Data Analysis 



In this section, we confirm the results in the previous section employing data on the assessed 
value of land in 1974-2007 Japan. In Japan, land is a very important asset which is distinguished 
from buildings. The assessed value of land indicates the standard land prices evaluated by 
Ministry of Land, Infrastructure and Transport Japan. The investigation is undertaken on each 
piece of land assessed once a year according to the posted land price system from 1970.^ 

The probability distribution functions of land prices are shown in Fig. 6. The number of 
data points is "14,570", "15,010", "15,010", "15,010", "15,580", "16,480", "17,030", "17,380" 
and "17,600" in 1974-1982, respectively. In 1983-1991, "16,975", "16,975", "16,975", "16,635", 
"16,635", "16,820", "16,840", "16,865" and "16,892", respectively In 1992-1999, "17,115", 
"20,555", "26,000", "30,000", "30,000", "30,300", "30,600" and "30,800", respectively. In 2000- 
2007, "31,000", "31,000", "31,520", "31,866", "31,866", "31,230", "31,230" and "30,000", re- 
spectively. 

In each figure, the power-law is observed in the large scale region. In addition, Pareto index 
H varies annually and changes significantly before and after bubble years (1986-1991). This is 
represented in Fig. 7 where each Pareto index /j, is estimated in the range of land prices from 
2 X 10^ to 10^ yen/m^. In Fig. 6, the log-normal distribution is also observed in each middle 
scale region. The standard deviation a varies annually and changes significantly in bubble years. 
This is also represented in Fig. 7 where each standard deviation a is estimated in the range from 
5 X 10^ to 3.17 X 10^ yen/m2. 

It is well known that the Pareto index for firm sizes hardly changes (Fig. 2). In such a 
database, the detailed balance (7) is observed [11] (Fig. 1). On the other hand, the Pareto index 
for the assessed value of land varies annually. This suggests that the system is not stable and the 
detailed balance does not hold. Actually most of the period, the detailed balance is not observed 
in the scatter plot of all pieces of land assessed in the database (Figs. 8 and 9 for instance).^ 

In the previous section, we have proposed the detailed quasi-balance (15) in an ideal qua- 
sistatic system. In each scatter plot of all pieces of land assessed in the database, the parameter 
9 is measured and the result is shown in Fig. 11. Here 6 is estimated in the following two regions. 
One is the large scale region between 2 x 10"'' and 10^ yen/m^ where Pareto index /i is estimated. 
The other is the middle scale one between 5 x 10'^ and 3.17 x 10^ yen/m^ where the standard 
deviation a of the log- normal distribution is estimated. In Fig. 11, we represent two parameters 
by and respectively. 

By using these parameters, we confirm not only the relation between the ratio of ^ + 1 and 

(27) in the large scale region but also the relation between the ratio of a and (30) in 
the middle scale one (Figs. 12 and 13). This warrants approximations assumed in the previous 

'^We exclude data in 1970-1973, the number of which is insufficient. 

*0f course, the detailed balance is observed approximately in the case that the system is almost stable and 
the Pareto index hardly varies (Fig. 10 for instance). 
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section. The derivation and the data analysis are consistent. 

5 Conclusion 

In this study, by using no model, we have derived a quasistatically varying log-normal dis- 
tribution from a Non-Gibrat's law under the detailed quasi-balance. In the derivation, we have 
employed two approximations. One is that the probability density function of the growth rate is 
described as tent-shaped exponential functions. The other is that the value of the origin of the 
growth rate distribution is constant. Even under the detailed quasi-balance, the first approxi- 
mation uniquely fixes a Non-Gibrat's law to be the same expression under the detailed balance. 
Together with the second approximation, the resultant distribution is described as power-law 
with varying Pareto index in the large scale region. In the middle scale region, the distribution is 
reduced to the quasistatic log-normal distribution with the varying standard deviation. Notice 
that not only the change of Pareto index n but also the change of the standard deviation a 
depends on a parameter 6 of the detailed quasi-balance. 

Employing empirical data on the assessed value of land in 1974-2007 Japan, we have con- 
firmed these analytic results. In the scatter plot of all pieces of land assessed in the database, the 
parameter 9 of the detailed quasi-balance is measured in the following two regions. One is the 
large scale region where Pareto index n is estimated. The other is the middle scale region where 
the standard deviation a is estimated. We have observed that two parameters On measured in 
the large scale region and measured in the middle scale one are in good agreement with the 
ratio oi fi + 1 and cr, respectively. 

Intriguingly, it is observed that the change of Pareto index influences the change of the 
standard deviation in the opposite direction (Fig. 7).^ In this database, the change of the 
distribution in the large scale region propagates in the middle scale one. This phenomenon is 
intelligible by two slopes and of the regression lines in the scatter plot (Fig. 11). 

In order to maJie these discussions more precisely, we need to investigate wealth of data in 
unstable state for a long period. By analyzing the database, we directly observe the Non-Gibrat's 
law which defines the large and middle scale regions. In addition, we should take a statistical 
test for the symmetry in the two arguments of Pi2{xi,X2) to confirm the detailed quasi-balance 
directly. These accurate data analyses are imperative to realize the applications of the study in 
this paper, such as credit risk management and so forth. 
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Figure 1: The scatter plot of all firms in the database, the profits of which in 2004 (xi) and 
2005 {x2) exceeded 0, xi > and X2 > 0. The number of firms is "232,497 ". 
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Figure 2: Probability density functions of profits in the database. In the large scale region, 
Pareto's law is observed. Each Pareto index is estimated to be nearly 1. 
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Figure 4: The relation between the lower bound 
of each bin xi and t±{xi) with respect to 
the profits growth rate from 2004 to 2005. 
Prom the left, each data point represents n = 
l,2,---,20. 



Figure 5: The probability distribution function 
of profits P{xi) for firms, the profits of which 
in 2004 (xi) and 2005 (X2) exceeded 0, xi > 
and X2 > 0. 
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Figure 6: Probability distributions of the assessed value of land in 1974-2007. In each figure, 
data points are equally spaced in logarithm of land price. 
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Figure 7: Annual change of Pareto index /x and the standard deviation a of the log-normal 
distribution from 1974 to 2007. 
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Figure 8: The scatter plot of all pieces of land assessed in the database, the values of which in 
1995 (xi) and 1996 (^2) exceeded 10^ yen/m^. The number of data points is "29,590". 




Figure 9: The scatter plot of all pieces of land assessed in the database, the values of which in 
1995 {xi) and 1996 (^2) exceeded 10^ yen/m^. The number of data points is "29,692". 
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Figure 10: The scatter plot of all pieces of land assessed in the database, the values of which in 
1978 (xi) and 1979 (0:2) exceeded 10^ yen/m^. The number of data points is "13,431". 
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Figure 12: Annual change of and (/xi + l)/(/X2 + 1) in the year {xi,X2) = (1974,1975) - 
(2006,2007). 




